Dispersive Lineshape Theory 
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Motivated by recent experiments we consider a stochastic lineshape theory for the case when 
the underlying process obeys power-law statistics, based on a generalized Anderson-Kubo oscillator 
model. We derive an analytical expression for the lineshape and find rich type of behaviors when 
compared with the standard theory, for example, new type of resonances and narrowing phenomena. 
We show that the lineshape is extremely sensitive to the way the system is prepared at time t — 
and discuss the problem of stationarity. 
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Since its first introduction by Anderson and Kubo(AK) 
lineshape theory based on a stochastic approach has 
had wide applications in condensed phase spectroscopy 
ranging from magnetic resonance spectroscopy |^,0| to 
the recently developed single molecule spectroscopy . 
The AK approach as well as other standard approaches 
(e. g. Bloch equation) is based on the Markovian assump- 
tion. These approaches have been applied with great suc- 
cess mainly to ensemble averaged measurements. How- 
ever, single molecule spectroscopy [HJ^ has revealed that 
in some cases the underlying dynamics in the condensed 
phase is highly non- Markovian [|-^. 

In recent fluorescence intcrmittency studies on single 
quantum dot systems [Ei 15| it was found that single 



quantum dots undergo transitions between bright and 
dark states during the measurement, and the sojourn 
times of both bright and dark states are distributed ac- 
cording to power-law, t~'^/^, in contrast to an exponen- 
tial as in the Poisson process. Here we are interested in 
the stochastic lineshape theory for such non-Markovian 
processes with emphasis on power-law behavior of the 
underlying dynamics. The introduction of a power-law 
process in the lineshape problem appears quite naturally. 
Consider the following model as an example. The static 
chromophore C at the origin interacts with the single 
perturber P via a short range interaction depending on 
r, the distance between C and P. C absorbs light with a 
frequency when r < and with a frequency oj'q ^ ujq 
when r > Vc- When P moves away from C it performs a 
random walk (assumed one dimensional for simplicity). 
Since the probability density of the first return time in 
the one dimensional random walk follows t""^/^ p6| , p7[ , 
the chromophore absorbs light at frequency lo'q with times 
distributed according to a power-law distribution. 

An important issue in a power-law stochastic process 
is stationarity, which is of concern due to a very broad 
temporal distribution for underlying processes. The sig- 
nificant difference between stationary and nonstationary 
cases has manifested itself in many physical problems, 
e. g. transport properties in disordered materials [ p^ and 
power-spectra in chaotic systems [|l9|. We take into ac- 



count the stationarity issue fully and show that the line- 
shape is a very sensitive measure of stationarity when the 
underlying dynamics obeys power-law statistics. 

The stochastic lineshape theory is based on the equa- 
tion of motion for the transition dipole, /i(t) = iLo{t)^{t), 
where u){t) is the stochastic frequency of the oscilla- 
tor. The dynamical quantity which determines the line- 
shape is the relaxation function <I>(i, ig) = (/^(^/^(^o)*), 
where the average is taken over all the possible real- 
izations of the underlying stochastic process and we 
have set (|Ai(io)P) — 1- From the equation of motion 
we can calculate the relaxation function as <!>(<, ip) — 
(^exp (^i J*^dTU}{T)j^. When the process is assumed 

to be stationary, ^{t,to) —^{t — to), then the line- 
shape I{w) can be calculated as the Fourier transform 
of the relaxation function by making use of the Wiener- 
Khintchine(WK) theorem p^. 
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where the symmetry of $(t), <&(— i) = ^*{t), has been 
used and the Laplace transform of z{t) denoted by z(s) — 
C{z{t)} and e ^ 0+. 

We assume that the underlying process is a renewal 
process as in the AK approach. To make the model as 
simple as possible, we consider a two state model. The 
transition frequency uj{t) of the chromophore can take the 
value of either — wp or +loo depending on the perturber 
state, |-|-) or |— ), respectively. Each alternating path be- 
tween the states |-|-) and |— ) of the perturber leads to 
a stochastic realization of chromophore frequency mod- 
ulation, and it is characterized by a sequence of sojourn 
times in the states |-|-) and |— ). The sojourn times in 
the states |±), t±, are assumed as mutually independent, 
identically distributed random variables described by the 
probability density functions (PDFs), h±{t±). The origi- 
nal AK process amounts to the exponential sojourn time 
PDF. We do not assume any specific functional forms for 
the sojourn time probability densities from the begin- 
ning, but are mainly interested in the process where the 
sojourn times are distributed with long time power-law 
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tails, (a > 0). 

Assuming stationarity, which is justified if the process 
has been going on for long times before the beginning of 
observation, the sojourn time PDFs for the first transi- 
tion event is given by f±{t±) [different than h±(t±)] by 
a standard argument pl[|. 
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where the mean sojourn time t±— dtth±{t) is as- 
sumed to be finite. When rj- — > cx) the concept of sta- 
tionarity breaks down. For the Poissonian case h±{t±) = 
exp{—t±/T±), so we have f± = h±. Therefore, sta- 
tionarity is naturally satisfied in the Poissonian process. 
However, the non-Poissonian process in which we are in- 
terested will not be stationary if we simply set f± ^ h±, 
and the WK theorem therefore does not hold. 

The conditional relaxation functions = -I-,—) 

are defined over the stochastic paths that start from the 
state \i) at time and end with the state \j) at time 
t. Here we give a sketch of the derivation for $++ by 
summing all the possible stochastic paths that start from 
and end at the state |-|-). Along a particular path if no 
transition is ever made until t, then the contribution of 
this path to is given by F^{t+)e~'''^°^+ in the time 
domain, where F±{t±) — J^^ dTf±{T) is the survival 
probability of the states |±) for the first event. This 
contribution will amount to -F+(s -|- ituo) in the Laplace 
domain. The next possible paths are those which make 
the first transition to the state |— ) at jump back to 
the state |-|-) after remaining at the state |— ) for time t2-, 
and stay at the state |-|-) until time t. The contribution 
of these to $++(f) is given by dti^ dt2-J^ dt'i,^ 
X /+ )e-^'^°*i+ h-{t2- )e"^°*^- H+ {t3+)e-"^°*'>+ with 
the constraint ti+ + t2- + t3+ = t, and H±{t±) — 
J^^ dTh±{T) —T±f±{t±) is the survival probability corre- 
sponding to h±{t±). In the Laplace domain this will read 
as + iujo)h^ (s — zciJo)-ff+(s + io^o) by the convolution 
theorem. Summing all the possible stochastic paths, we 
have 



^++{s)^F+ + f+h^{l + h+h^ 
f+h-H+ 
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where x± = x±{s ± ioJo) with x being f,h,F, or H. In a 
similar way, we have 



<f+_(s) 



1 - h+h^ 



(4) 



and $__(s) and <t |-(s) are obtained from <I>-|__(-(s) and 

$^ (s) by exchanging -I- and — . 

The total relaxation function can be calculated from 
the conditional relaxation functions. 
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with the initial distribution of the perturber state given 
by P± — T±/J2i=± '^i- Then from Eq. (|l]) the lineshape 
is given by 



= -Re 
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where z± ^ iuj ± icuQ and h± ~ h±{z±), and we have ex- 
pressed f±,F±, and H± in terms of h±. This is the final 
expression of the lineshape function for the stochastic os- 
cillator undergoing the two state frequency modulation. 

It is our aim here to show that the theory exhibits a 
very strong sensitivity on the choice of PDF for the first 
event. This becomes important for experimental situa- 
tions when it is not always clear if the underlying process 
is stationary or not. For this purpose we define a quasi 
lineshape Ins{'^) by replacing f± with h± in the deriva- 
tion of Eq. (0), 
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Note that for Poissonian case I^si'^) = a-nd for or- 
dinary processes (i.e. when all moments of h±{t) exist) 
Ipfsi^) ~ However, as we show here, strong sen- 

sitivity on the first event is recovered for the dispersive 
case. We emphasize that Eq. (^) is not a lineshape be- 
cause the underlying process is not stationary. In cer- 
tain physical situations, however, the underlying pro- 
cess is not an ongoing process, but has been initiated 
by a measurement itself, for instance, the blinking pro- 
cess in the quantum dot experiment mentioned before. 
In this case one can calculate the relaxation function 
^Ns{t) = ^{t, 0) for the nonstationary stochastic process 
and /jvs(a;) is obtained as the complex Laplace trans- 
form of a well-defined nonstationary relaxation function. 

The original AK model is recovered from Eq. (|^) by 
choosing an exponential sojourn time PDF. 0J^. We 
first consider the sojourn time PDFs which have finite 
first moments, t± < oo, but divergent second moments. 
As a representative of this class, we use the following 
form, 



h±{t±) = (r3/27r4)'/%xp(-r±/2t±). 



(8) 
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In this case, h±(t) decays as t at long times, thus the 
first moment exists, but the second moment diverges. 
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FIG. 1. The lineshapes calculated from /i3/2(t) for the 
stationary and the nonstationary cases are compared with the 
AK case from slow to fast modulation limits. The resonance 
frequency is chosen as tJo = 1- The correlation time is chosen 
as uJoT — 20, 4, 0.4, and 0.08 in (a)-(d), respectively. 

In Fig. we have compared the stationary and non- 
stationary lineshapes for h±(t) in Eq. (||) and for the 
AK case. We have chosen /13/2 = h+ = h- by setting 
r+ = T_ = r. The correlation time of the perturber dy- 
namics is varied from slow {loot ^ 1) to fast {loqt <C 1) 
modulation cases in (a)-(d). For the AK case, the well 
known phenomenon of motional narrowing is shown: in 
the slow modulation case we see two peaks at a; = iwo 
while in the fast modulation case we observe a single peak 
at w = 0. For the case the stationary and nonsta- 
tionary cases show very different behaviors as the correla- 
tion time is decreased; thus, the first event in the underly- 
ing random process has a strong effect on the lineshapes. 
In addition, new phenomena are found for the station- 
ary lineshape in the fast modulation cases (Fig. |l](c) and 
(d)): three distinct peaks are observed for the stationary 
case. 

The new peaks we observe in Fig. |l|(c) and (d) at 
LD — ibwo for the stationary lineshape result from the first 
event in the stochastic process w(t). The probability for 
the perturber remaining at the initial state is governed 
by the long time tail in the sojourn time PDF. Due to the 
stationarity condition in Eq. (^) the survival probability 
for the first event decays more slowly for the stationary 
case(~ t^~°') than for the nonstationary case(~ t~°'), 
where a = 3/2 in this example. Therefore, the stationary 
case effectively requires the perturber to remain at the 
initial state until much longer times than the nonstation- 



arity case, resulting in the enhanced peaks at w = iwo. 
This is why we observe new peaks not present in the 
standard Poissonian case. 

As the next example we consider the one-sided Levy 
density as the sojourn time PDF jl7|, h{s) = La{s/r) — 
exp(— (s/r)"), with < a < 1, and r being a coefficient 
with an inverse time dimension. It is well known that 
the Levy PDF decays algebraically at long times ^ 1, 
La{t) ~ and thus all the moments of La{rt) in- 

cluding the first moment diverge Therefore, there 

is no microscopic timescale for this PDF, and the form 
of f±{t) given in Eq. (^) cannot be applied. However, in 
realistic situations, it is unlikely to have power-law be- 
havior for an infinitely long time, but rather, it is likely 
to have a finite cut-off time provided by, for example, 
the lifetime of a molecule. Therefore, it is natural to in- 
troduce a cut-off time tc such that the algebraic decay 
is valid during time interval <^ t <^ tc- We intro- 
duce an exponential cut-off function for the convenience 
of an analytical treatment. Now the sojourn time PDF 
is given by h±{t) = N±e~*-^*"' La{r±t), where Af± is the 
proper normalization constant depending on the cut-off 
time. Then the Laplace domain expression of h±{t) can 
be written as 

h±{s) = exp [(r±t,)-" {1 - (1 + st,)"}] . (9) 

Then f±{t) is given from Eq. (H) with the mean be- 
ing T± = atc/{r±tc)°' . Note that in the limit tc ^ 00 
the Levy PDF without cut-off is recovered as h±{s) — 
exp(— (s/r±)") and t± diverges. 
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FIG. 2. The lineshapes are calculated with the Levy PDF 
cases with two different values of the Levy index, a. The 
other parameters are chosen as tJo = 1 and r± — 1, except 
for the cut-off time which is varied, (a) a — 0.3, stationary, 
(b) a = 0.3, nonstationary, (c) a — 0.8, stationary, and (d) 
a = 0.8, nonstationary cases. 
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In Fig. 1^ we have investigated the effect of the cut- 
off time in the Levy PDF case both in the stationary 
and the nonstationary cases. When a = 0.3 (Fig. | 
(a) and (b)), both the stationary and the nonstation- 
ary Levy hneshapes show distinct peaks at w = ±wo. 
As the cut-off time is increased, there is httle change in 
both hneshapes other than becoming narrower. When 
a = 0.8 (Fig. H (c) and (d)), there appears a new peak 
near lo = Q 'm addition to the two resonance peaks. This 
is a new type of the narrowing behavior which is absent in 
the Poissonian case in that it is controlled by the power- 
law index a rather than the correlation time r (as in 
the Poissonian case), and is termed power-law narrowing. 
Also, as the cut-off time is increased, the central peak in 
the lineshape for the stationary case diminishes while it 
remains in the nonstationary case. This is because in 
the stationary case, as the cut-off time is increased the 
first event will dominate the probability weight in the 
stochastic paths of the perturber dynamics. The differ- 
ence between the stationary and nonstationary cases is 
therefore more significant in the Levy case than in the 
case /i3/2(t) given in Fig. 1. 

To investigate the power-law narrowing we consider 
the limit tc ^ oo. In this limit the stationary lineshape 
approaches two delta functions, = P+5UjJ + ujq) + 

P-^5{bj — bJi)), since the second term in the Eq. (|l|) vanishes 
as T± —t oo. The nonstationary case, however, yields in 
the limit r± ^ 



lim iNsi^^) 

tc — *00 
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r]x'^ + (77a;") -1 + 2cos(7ra) ' 
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for |ti;| < loq and vanishes when \lu\ > luq. This ex- 
pression has been obtained from Eq. (^ by taking the 
small frequency limit of the sojourn time PDF, h±{s) = 
1 - (s/r±)" -I- Note that Eq. is not limited 

to Levy PDF, but valid for any PDF with ^"(1+") tail 
(0 < a < 1). Here, x = {luq + uj) / {ojq — uj) is a dimension- 
less frequency, and 77 = limt^^oo(p+/p-) = (^-/^+)" 
is the asymmetry parameter. This function shows a 
very asymmetric power law singularities at a; = ii^a, 
{ujQ ± uj)^^°' depending on a. It is worthwhile to men- 
tion that such a strong asymmetric lineshape has been 
encountered in the problem of the X-ray edge absorption 
of metals In the symmetric case (77 = 1) Eq. ( p^ ) 

reduces to a simpler expression, and in this case there 
exists a critical value of a below which the lineshape is 
concave and above which convex at cj = 0, which is given 
by ac = cos(7rac/2) — 0.5946 • • •. 

In Fig. ^ we have confirmed this finding by plotting 
the nonstationary hneshapes for the Levy PDF case with 
rj = 1 and finite, large r± and t^.- The nonstationary case 
in Fig. H (b) shows the concave-to-convex transition at 
the critical value of a as predicted. For the stationary 
lineshape, similar kind of behavior can be observed in 
Fig. H (a), however, Uc now depends on tc- 
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FIG. 3. The power-law narrowing for the lineshape with 
Levy PDF is shown as the Levy index a is changed. The 
other parameters are chosen as r± — 100 and tc = 10*. (a) 
stationary and (b) nonstationary cases. 



In summary, we have generalized the stochastic line- 
shape theory to arbitrary renewal processes. Compared 
with the standard theory, we have found a variety of new 
phenomena in the hneshapes, such as new peaks and nar- 
rowing behaviors. The issue of the stationarity has been 
considered and the strong sensitivity of the lineshape to 
the first event in the stochastic trajectory was found. One 
of many extensions of this work is to consider the ran- 
dom frequency modulation of stochastic oscillator among 
N > 2 values, and will be published elsewhere. This work 
was supported by NSF. 
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